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ABSTRACT 

We report optical observations of the Luminous Blue Variable (LBV) HR Carinae which show that 
the star has reached a visual minimum phase in 2009. More importantly, we detected absorptions 
due to Si IV AA4088-4116. To match their observed line profiles from 2009 May, a high rotational 
velocity of Wrot — 150 ± 20 kms"'^ is needed (assuming an inclination angle of 30°), implying that HR 
Car rotates at ~ 0.88 ± 0.2 of its critical velocity for break-up (fcrit)- Our results suggest that fast 
rotation is typical in all strong- variable, bona-fide galactic LBVs, which present S Dor-type variability. 
Strong-variable LBVs are located in a well-defined region of the HR diagram during visual minimum 
(the "LBV minimum instability strip"). We suggest this region corresponds to where Vcrit is reached. 
To the left of this strip, a forbidden zone with Wrot/wcrit > 1 is present, explaining why no LBVs are 
detected in this zone. Since dormant/ex LBVs like P Cygni and HD 168625 have low Urot, we propose 
that LBVs can be separated in two groups: fast-rotating, strong- variable stars showing S-Dor cycles 
(such as AG Car and HR Car) and slow-rotating stars with much less variability (such as P Cygni and 
HD 168625). We speculate that SN progenitors which had S-Dor cycles before exploding (such as in 
SN 2001ig, SN 2003bg, and SN 2005gj) could have been fast rotators. We suggest that the potential 
difficulty of fast-rotating Galactic LBVs to lose angular momentum is an additional evidence that 
such stars could explode during the LBV phase. 

Subject headings: stars: atmospheres — stars: mass loss — stars: variables: other — supergiants — 
stars: individual (HR Carinae) — stars: rotation 



1. INTRODUCTION 

Massive stars are the main contributors to the in- 
put of ionizing photons, energy, and momentum into 
the interstellar medium and are responsible for a signif- 
icant fraction of the chemical enrichment of their host 
galaxy. Evolutionary models predict the existence of 
a short-lived, transitional stage, usually re ferred to as 
the Ltiminous Blue Variable (LBV) phase (|Contilfl98l 
iHumphrevs fc Davidsonlll994f ). characterized by a high 

mass-loss rate (M ~ 10~^ - 10"^ Moyr"^). In the 
current picture of stellar evolution, LBVs are rapidly- 
evolving massive stars in the transitory phase from be- 
ing an 0-type star burning hydrogen in its core to be- 
coming a Wolf-Rayet (WR) type, helium core-burning 
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star (Hum phrevs fc Davidson' '1994'; iMaeder fc MevneB 
20Q0b: M evnet fc Maeder 2000 . 2003,). 

There is strong evidence that very massive stars lose 
the bulk of thei r mass through giant o utbursts during 
the LBV phase (|Smith fc Owockill2006h . similar to the 
one that occurred in Eta Carinae in the 1840's. Surpris- 
ingly, recent works suggest that some core -colla pse super 
novae (SNe) have LBV progenit or s fe.g.jKota k fc Vin 
20061: lSmith "erar 1 20071 [2008allbl: iGal -Yam e t al.. ,200 ' 



Trundle et al. 2008; iGal-Yam fc Leon ard 2003), which 



not only dramatically enhances the cosmological impor- 
tance of LBVs but also poses a great challenge to the 
current paradigm of massive star evolution. 

LBVs are rare: only a dozen or so are known in 
the Galaxy, and roughly only three of them, namely 
AG Carinae, HR Carinae, a nd WRA 751, are con- 
firmed strong-variable LBVs (Humphrevs fc DavidsonI 
IT99I Ivan GenderenI l200lt iClark et al.l [2OO5I1 These 
stars have been observed for a sufficiently long time to 
establish a record of strong photometric and spectro- 
scopic variability which shows, undoubtedly, that both 
strong S Dor-type variability and a circumstellar nebula 
are present. The spectroscopic and photometric behav- 
ior of these strong-variable LBVs are markedly differ- 
ent compared to other stars which are usually regarded 
as dormant LBVs, such as P Cygni and Eta Carinae 
(Ivan GenderenI [20?m . 

AG Car, HR Car, and WRA 751 are unstable and 
present cyclical strong S Dor-type variability character- 
ized by irregular visual magnitude changes on timescales 
of decades, with a typical amplitude of ISV ~ 1 — 2 mag, 
and corresponding changes in effe ctive t emperature and 
hydrostatic radius (jvan GenderenI [20011 ). During visual 
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Fig. 1.— ASAS-3 CCD V-band lightcurvo of HR Carinae from 
2001-2009. The arrows indicate epochs when spectroscopic data 
were available during minimum. 

minimum, the star is typically hot, while at visual max- 
imum, a cooler effe ctive temperature is obtained (e.g., 
Ivan Gendered [r982l ) . How the S Dor- type variability 
relates to the powerful giant eruptions is not clear, al- 
though it could be possible that a relatively large amount 
of stellar mass, which is not ejected fro m the star, is tak- 
ing part in the S Dor- type variability (jGroh et al.ll2009l : 
hereafter GHD09). This would suggest t hat the S Dor- 
type variability is a failed giant eruption ([GHDOgf ) . 

At least for AG Car, two additional pieces of the puz- 
zle are known: a significant reduction (~ 50%) in the 
inferred bolometric luminosity fr om visua l minimum to 
maximum has been determined (jGHDOQl ). and a high 
rotational velocit y has been obtained during minimum 
(|Groh et al.ll2006l : hereafter GHD06). Whether the same 
is true for other bona-fide LBVs, such as HR Car, still 
remains to be seen. 

In this Letter we report photometric and spectroscopic 
observations of the bona-fide LBV HR Car obtained dur- 
ing 2009 (§[2|), showing that the star has entered a visual 
minimum phase of the S-Dor cycle. More importantly, we 
detected absorptions related to the high-excitation lines 
of Si IV AA4088-4116 in HR Car. The resuhs of our non- 
LTE spectroscopic modeling using the radiative trans- 
fer code CMFGEN are presented in § [3] They suggest 
a photospheric nature for these lines, that HR Car has 
high rotational velocity (wrot), and that Ui-ot is very high 
compared to the critical rotational velocity for break- 
up (^^crit)- In § m we discuss the importance of such a 
detection for understanding the LBV phase and the sub- 
sequent evolution of very massive stars. 

2. OBSERVATIONS 

Optical CCD F-band photometr y of HR Car was com- 
piled from the ASAS-3 database froim anskill2002f ) and 
averaged in bins of 15 days, with only high-quality (grade 
A) data being used. The ASAS-3 magnitudes were mea- 
sured using a large aperture (75"), but HR Car is by far 
the dominant source of light in the field. We estimate 
that the typical error of the 15-day averaged y-band 
magnitude is 0.03 mag. Figure [1] presents the F-band 
lightcurve of HR Car from 2001-2009, showing that a 
significant increase of 1.5 mag occurred in this time span 
and that HR Car has entered a visual minimum phase of 
the S-Dor cycle. 

Optical spectra of HR Car around the Si IV lines 
at AA4088 — 4116 were obtained using several facilities 
in the southern hemisphere. Medium resolution (i? 2± 
12, 000) spectroscopic observations were made at the 1.6- 
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Fig. 2.— Optical spectra of HR Car around Si IV AA4088-4116. 
From bottom to top, we show spectra obtained from 2009 January 
13 until 2009 May 27. The blue wing of Si IV A4088 was hit by 
cosmic rays on 2009 January 18th and is shown by the dashed line. 

m telescope of the Observatorio Pico dos Dias (OPD- 
LNA, Brazil) and at the 2.15-m telescope of the Com- 
plejo Astronomico El Leoncito (CASLEO, Argentina). 
High-resolution spectra (i? ~ 48, 000) were gathered us- 
ing the Fiber-fed Extended Range Optical Spectrograph 
(FEROS; Kaufer et al.l ll999) mounted on the 2.2-m tele- 
scope of the Max-Planck-Gesellschaft/European South- 
ern Observatory (MPG/ESO, Chile). The reduction of 
the spectroscopic data, including bias and sky subtrac- 
tion, flat-ficlding, spectrum extraction, wavelength cali- 
bration, correction to the heliocentric frame, and contin- 
uum normalization, was carried out using sta ndard IRAF 
routi nes. A systemic velocity of —10 kms""'^ (jWeis et al.l 
[1991 was subtracted from the observed spectra. Table [1] 
summarizes the observational data of HR Car available 
during visual minimum. 

Figure [2] presents a montage of the observed Si IV 
AA4088 - 4116 line profiles of HR Car during our 2009 
monitoring. The use of different instruments over 5 
months confirms that the detection is not transient and, 
undoubtedly, absorptions due to Si IV AA4088 — 4116 
lines are typical in HR Car during visual minimum. 

3. MODELING AND RESULTS 

In the case of the prototypical LBV AG Car, the rela- 
tively broad absorptions due to Si IV AA4088 — 4116 are 
caused b y the extremely fast rotation, Urot/'i'crit ^ 0.86 
(|GHD06D . Here we also use these Si IV lines to derive 
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TABLE 1 

Journal of optical spectroscopy of HR Car, during the 2009 visual minimum 



Date V Telescope/Instrument Spectral Range (A) R 



2009 January 13 8.74 2.15-m CASLEO/REOSC 3650-6130 12,000 

2009 January 18-19-23 8.77 2.15-m CASLEO/REOSC 3650-6130 12,000 

2009 April 18 8.66 1.6-m OPD-LNA/Coude 3880-4390 12,000 

2009 May 03-27 8.67 2.2-m MPG-ESO/FEROS 3600-9200 48,000 



the rotational velocity of HR Car , perform ing a detailed 
spectroscopic analysis similar to GHD06, using the ra- 
diative transfer code CMFGEN (HiUier fc Miller 199a) 
and th e two-dimensional (2-D) code of iBusche fc Hillien 
((20051 hereafter BH05). 

CMFGEN solves the statistical equilibrium equations 
simultaneously to the radiative transfer in the co-moving 
frame of the outflow, in spherical symmetry, and assum- 
ing steady state. The code computes line and contin- 
uum formation in the non-LTE regime, and each model 
is specified by the stellar luminosity L^, and effective tem- 
perature Toff®, the wind mass-loss rate M, volume fill- 
ing factor /, and terminal velocity Woo, and the chem- 
ical abundances Zi of the included species. Since the 
momentum equation of the wind is not solved, a beta- 
type velocity law is adopted and modified at depth to 
smoothly match a hydrostatic structure at R^,. The hy- 
drostatic equation is iterated in order to produce a quasi- 
hydrostatic structure which extends inwards until the in- 
ner boundary of the model (at a Rosseland optical depth 
of 100). Line blanketing is fully taken into account, and 
our models include tens of thousands of spectral lines in 
non-LTE. The atomic mode l used fo r HR Car is similar 
to the one used for AG Car (|GHD09f ) and included lines 
of H, He, C, N, O, Na, Mg, Al, Si, Cr, Mn, Fe, Co, and 
Ni. 

We assumed a depth-independent turbulent velocity 
of 20kms~^ in the computation of the observed spec- 
trum, which is an upper limit, since a higher turbulent 
velocity would yield a detectable redshift, which is not 
observed, of the emission lines with high opacity, such 
as the stronger lines of H I and He I. Therefore, mi- 
croturbulence cannot explain the broadening seen in the 
high-ionization Si IV AA 4088-4116 lines, and we suggest 
these lines to be rotationally broadened (note that other 
broadening mechanisms, such as convection or other kind 
of surface activity, have yet to be detected in LBVs and, 
thus, are not expected to be present in AG Car, although 
our data do not rule them out). 

The effects of rot ation o n the line profiles were in- 
vestigated using the iBHOSi code, which allows the com- 
putation of the spectrum in 2-D geometry for a given 
inclination angle i. No equatorial or polar density en- 
hancements were taken into account, since our goal was 
to investigate the broadening of the photospheric lines. 
Therefore, we computed a grid of synthetic line pro- 
files with a stellar rotational velocity ranging from to 
200 kms~^ in increments of 5 kms""'^, and a best fit was 
obtained by visual inspection. We assume that the stellar 
rotational axis of HR Car is aligned with the polar axis 
of its bipolar nebula, which is inclined at i = 30° ±10° 

* The effective temperature is defined as the temperature where 
the Rosseland optical depth is 2/3, i.e., T^g = T(-7-ross = 2/3). 



(|Nota et al.lll997[ ). 

In a dense stellar wind, ionization stratification oc- 
curs, and lines from ions of different ionization poten- 
tials will form at different distances r from the star. The 
Si IV AA4088 — 4116 lines are the highest ionization lines 
present in the optical spectrum of HR Car and, therefore, 
should be formed closest to the stellar photosphere. Since 
the azimuthal velocity component is proportional to r~^, 
Si IV AA4088 — 4116 are ideal for deriving the rotational 
velocity of HR Car. Our modeling shows that other opti- 
cal lines which are ocasionally seen in absorption in LBV 
spectra, such as N II AA4601 - 4607 - 4630 - 4643 and 
Si III AA4553 - 4558 - 4575, are actually formed in the 
wind of HR Car, far from the stellar photosphere. As a 
consequence, their broadening is dorn inat ed by the wind 
velocity field (see discussion in [BH05l and lGHD06h . 

Figure [3] compares the absorption profiles of Si IV 
AA 4088-4116 observed in 2009 May 03 (which has the 
highest signal-to-noise ratio and spectral resolution of 
our sample) with our best CMFGEN model, which has 
the following parameters: = 5. x 10^ Lp, (assum- 
ing Ay = 3.1 and d = 5 kpc. Ivan Genderen et al.l 
IT99lt iHutsemekers fc van Dro3 |199lD, R* = 70 Rq, 
i?phot = 75 Rq, Toff = 17900 K, He/H=0.4 (by number), 

Voo = 120 kms-\ /3 = 1, M = 7 X lO'^ Mqji'^ and 
/ = 0.25. In a forthcoming full paper we will present the 
results of the detailed spectroscopic analysis of HR Car 
during minimum, including diagnostics for each of the 
parameters assumed above. Nevertheless, for reasonable 
changes in the above parameters, the derived rotational 
velocity of HR Car changes only slightly and the conclu- 
sions of this Letter remain valid. 

Similar to AG Car (,GHD06i) . the observed absorption 
profiles of both Si IV lines are significantly broader and 
shallower than the non-rotating model. Our model pre- 
dicts that the Si IV lines are formed very close to the 
stellar photosphere, where the wind is still accelerating 
and, therefore, the line broadening caused by the wind 
velocity field is negligible. Consequently, most of the line 
broadening should be due to rotatio n. Acc ording to our 
rotating models computed with the IBH05I code, a rota- 
tional velocity of Wfot — 150 ± 20 kms~^ is required to 
reproduce the observed shape of the Si IV AA4088~4116 
fines in HR Car (Fig. ^. 

How does the rotational velocity of HR Car compare to 
its critical velocity for break -up? An estimation of fcrit 
can be calculated following iMaeder fc Mevnetl (|2000aD 
once the stellar mass M, Li,, iJ^, and the Eddington pa- 
rameter r are known. The critical parameter for ob- 
taining a precise value of ■yrot/wcrit is M, which can 
be estimated based on the evolutionary models from 
iMevnet fc Maedeil (|2003f ). Such models predict that, 
based on our derived value of ~ 5.0 x 10^ Lq, the 
initial mass of HR Car should have been roughly Mini — 
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Fig. 3. — Comparison between the Si IV AA 4088-4116 line pro- 
files of HR Car observed on 2009 May 3 rd with synthetic model 
spectra computed with the 2-D code from lBHOSI for different rota- 
tional velocities, assuming i = 30°: kms-i (green), 75 kms'^ 
(red), 150 kms"-*^ (blue), and 200 kms"-*^ (purple). 



50 ± 10 Mq. While there are no evolutionary models 
available in the literature specifically for Mini = 50 Mq, 
we estimate that stars such as HR Car will reach the 
LBV phase with M ~ 25 Mq, since evolutionary models 
predict that a 40 Mq star reaches the LBV phase with 
20 Mq, while a 60 Mq star reach es the LBV phase with 
30 Mq (|Mevnet fc MaedeHl2003D . 

Assuming an admittedly uncertain M — 25 Mq, our 
CMFGEN model of HR Car predicts F ~ 0.80, and we 
derive Wcrit — 170 ±20 kms~^ and_Urot/wcrit — 0.88 ±0.2. 
Therefore, similar to AG Car ()GHD06l ). HR Car is also 
likely rotating very close to, if not at, its critical velocity 
for break-up. 

4. DISCUSSION: STRONG-ACTIVE, BONA-FIDE 
GALACTIC LBVS ARE FAST ROTATORS 

Several theoretical works have long predicted that 
LBVs should have hi g h rot ati onal velocitie s (e.g . 
Humphrcvs & Davidson' '1994'; 'Langer et al.' '1999'; 
Maeder fc Mcvnet 2()00a; .Maeder fc Desiacques 2001; 
Dwarkadas fc Owockin2002t lAerts et all l200l ). but 
direct observational detection through the broadening 
of spectral lines has been elusive until recently, in 
particular because high-resolution observations of LBVs 
in the blue part of the optical spectrum are scarcely 
available during visual minimum. The detection of high 
rotational velocity in two of th e most variable LBVs in 
the Galaxy, AG Car (IGHD06I ) and now HR Car (§ E]), 
raises the obvious question: are all LBVs fast rotators? 

The presence of fast rotation in two out of three^ of 
the galactic bona-fide LBVs that have a circumstellar 
nebula and strong S Dor-type variability is probably not 
a coincidence, and we suggest that fast rotation is typical 
in such stars. Perhaps, fast rotation might even be a 
key ingredient for the instability that causes the S-Dor 
cycles in LBVs. If that is true, we should expect that 
fast rotation is also present during visual minimum in 
extragalactic LBVs, such as R127 and S Doradus in the 
LMC. An obvious test case for such a hypothesis would 
be R127, which h as been undergoing a visual minimum 
phase since 2007 (jWalborn et al.ll2008D . 

Figure [4] presents the position of AG Car and HR Car 

^ Unfortunately, no high-resolution observations of WRA 751 
during minimum are available to us to check if it is a fast rotator. 
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Fig. 4. — HR diagram showing the position of the strong- 
variable LBVs AG Car and HR Car during visual minimum, ac- 
cording to the updated stellar parameters determined in GHDQ^ 
and in this work. A revised position for the LBV minimum in- 
stability strip is provided (red dashed line) , and is considerably 
steeper than previous determinations l|Wolll 119891 : Ivan GendcrenI 
2001; Clark et al. 2005). Notice that AG Car and HR Car are 
both fast rotators during minimum, which suggests that the LBV 
minimum instability strip corresponds to the location where crit- 
ical rotation is reached for stron g-variable LBVs. The location of 
the Humphreys-Davidson limit IHumphrevs fc Davidsonlll994l) is 
shown (black dashed line). 



in the HR diagram basejl on their updated stellar pa- 
rameters obtained in , GHD09i and in this work. Pre- 
vious works have recognized that, during visual mini- 
mum, LBVs are located in a well-defined region of the 
HR di agram, the so-called "LBV minimum instability 
strip" (|Wolj|1989tlClark et al.ll2005D . Our results suggest 
that the LBV minimum instability strip is considerably 
steeper than previously determined, being characterized 
by log(L/Lo) = 4.54 log (Teff /K) - 13.61 (red dashed 
line in Fig. |4]). More importantly, we suggest that the 
LBV minimum instability strip corresponds to the region 
where critical rotation is reached for LB Vs with strong S- 
Dor type variability. When LBVs are evolving towards 
maximum, the star moves far from the LBV minimum 
instability strip (to the right in the HR d iagram), and 
Wrot/wcrit dccrcascs considerably (|GHD06l 1. The region 
in the HR diagram on the left side of the LBV minimum 
instability strip would be populated by unstable LBV 
stars with Vrot/i'crit > 1, which would make it a "for- 
bidden region" for LBVs. Indeed, no confirmed strong- 
variable LBV is s een in this region (jvan Gendered [2OOII : 
iClark et al.ll2005h . 

A corollary of the relationship between S Dor-type in- 
stability and fast rotation that we suggest in this Letter 
regards SNe which likely had an LBV progenitor. Some 
of these LBVs showed evidence for S Dor-type variability 
before exploding as a SN, such as the p rogenitors of SN 
2001ig, SN 20q3bg dKotak fc Vinkl l2006l). and SN 2005gj 
(|Trundle et al.ll2008D . It could be very likely that these 
massive stars were fast rotators before exploding. 
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What about dormant LBVs which have not presented 
evidence for strong S Dor-type variabiHty for a century or 
more? From these, only Eta Car has shown indirect evi- 
dence of fast rota tion, through the presence of a latitude - 
dependent wind (|Smith et al.l l2003: Wcig elt et al.l[2007f ). 
Interestingly, at least two other confirmed LBVs that do 
not show S Dor-type variability appear to rotate slowly. 
P Cygni has evidence of a relatively low rot ational ve- 
locity of 50 - 60 kms-i (jNaiarro et al.lll997D ^°. a prop- 
erty also shared by HD 168625, for which we measured 
tisini ~ 40 kms"^ based on optical photospheric N II 
lines and a similar methodology as in the present work 
(J. H. Groh 2009, in preparation). Assuming i ~ 60° 
(lSmithl[200l . we obtain Wrot — 45km s ^. For both P 
Cygni and HD 168625, the critical velocity for break-up 
should be of the order of 150-200 kms^-'^, meaning that 
^^rot/^^crit should be ~ 0.2-0.3 for these stars. 

We suggest that LBVs, at least in our Galaxy, can be 
separated into two groups: group I, comprised of fast- 
rotating, strong-variable LBVs such as AG Car and HR 
Car, and group II, composed of slow-rotating, weak vari- 
able LBVs such as P Cygni and HD 168625. While an 
LBV can in principle migrate from one subgroup to the 
other, for instance after a giant eruption, these LBV 
groups might also follow different evolutionary paths dur- 
ing and after the LBV phase, possibly yielding signifi- 
cantly different SN progenitors. 

In particular, fast-rotating LBVs might struggle to 
lose enough angular momentum through winds and gi- 
ant eruptions. Fast-rotating stars are supposed to have a 
latitu de-dependent, polar-enhanced wind (|Owocki et al. 
I1996D which, depending on the mass-loss rate, will not 
remove as much angular momentu m as in the ca se of 
a spherical wind (see discussion in lMaeded[2f)0l . In 
addition, LBVs might not lose enough angular mo- 
mentum during giant outbursts, similar to the one 
that occurred in Eta Car in the 1840's, since almost 

Assuming the rotational axis o f P Cygni is aligned with part of 
its nebula, which is seen at i ~ 45° l lSmith fc Hartiganll2003 ') . Note 
that this assumption is not crucial for the results of this paper. 

IMassevI 1)19801 ) suggested a high rotational velocity of ~ 
500 kms-l for the WN5 star HD 193077, based on broad absorp- 
tions due to high Balmer and He I lines. However, lAnnuM 111 9901 *) 



75% of the m ass was lost at latitudes higher than 45° 
(jSmithl l2006l ). Although current evolutionary mod- 
els predict that fast-rota ting LBVs should spin-down 
(jMevnet k Maedeil |2003[ ). the models do not include 
angular momentum and mass losses through a giant 
outburst, which i s a ke y factor during the LBV phase 
(jSmith fc Owockill2006D . In the Galaxy, there are at 
least three examples of stars which had a giant eruption, 
lost significantly amounts of mass (se veral M(^). but are 
still currently fast rotators: Eta Car ( Smith et al."2003|; 
iWeieelt et aLl[200l . AG Car (iGroh et al., 2006). and HR 
Car (this work). 

As a consequence, if fast-rotating LBVs are to evolve 
to WR stars without losing enough angular momen- 
tum, they would produce a population of significantly 
fast-rotating WR stars, which is yet unseen in the 
Galaxy^ ^ and not pred icted by evolutionary models 
(jMevnet fc Maedeil I2003D . An alternative evolutionary 
path would be that LBVs could fail to lose an additional 
significant amount of mass and thus, instead of becom- 
ing a Wolf-Rayet star, would dramatically explode as a 
core- collapse supernova anytime dur ing the LBV phase 
(e.g. ISmith et all 120071 : ISmithI l2008l ). We suggest that 
the potential difficulty of fast-rotating Galactic LBVs to 
lose angular momentum is an additional evidence that 
LBVs such as AG Car and HR Car could explode during 
the LBV phase. 
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proposed that HD193077 is a WR+O binary and that the broad 
absorption lines are formed in the O star. Even if the absorption 
lines were formed in the WR wind, they should be formed suffi- 
ciently far from the hydrostatic core that the broadening would be 
dominated by the wind acceleration. 
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